Nonlinear Acceleration of Sequential Fully Implicit (SFI) Method for Coupled Flow and Transport in Porous Media

The sequential fully implicit (SFI) method was introduced along with the development of the multiscale finite volume (MSFV) framework, and has received considerable attention in recent years. Each time step for SFI consists of an outer loop to solve the coupled system, in which there is one inner Newton loop to implicitly solve the pressure equation and another loop to implicitly solve the transport equations. Limited research has been conducted that deals with the outer coupling level to investigate the convergence performance. In this paper we extend the basic SFI method with several nonlinear acceleration techniques for improving the outer-loop convergence. Specifically, we consider numerical relaxation, quasi-Newton (QN) and Anderson acceleration (AA) methods. The acceleration techniques are adapted and studied for the first time within the context of SFI for coupled flow and transport in porous media. We reveal that the iterative form of SFI is equivalent to a nonlinear block Gauss-Seidel (BGS) process. The effectiveness of the acceleration techniques is demonstrated using several challenging examples. The results show that the basic SFI method is quite inefficient, suffering from slow convergence or even convergence failure. In order to better understand the behaviors of SFI, we carry out detailed analysis on the coupling mechanisms between the sub-problems. Compared with the basic SFI method, superior convergence performance is achieved by the acceleration techniques, which can resolve the convergence difficulties associated with various types of coupling effects. We show across a wide range of flow conditions that the acceleration techniques can stabilize the iterative process, and largely reduce the outer iteration count

MeshfreeFlowNet: A Physics-Constrained Deep Continuous Space-Time Super-Resolution Framework

We propose MeshfreeFlowNet, a novel deep learning-based super-resolution framework to generate continuous (grid-free) spatio-temporal solutions from the low-resolution inputs. While being computationally efficient, MeshfreeFlowNet accurately recovers the fine-scale quantities of interest. MeshfreeFlowNet allows for: (i) the output to be sampled at all spatio-temporal resolutions, (ii) a set of Partial Differential Equation (PDE) constraints to be imposed, and (iii) training on fixed-size inputs on arbitrarily sized spatio-temporal domains owing to its fully convolutional encoder. We empirically study the performance of MeshfreeFlowNet on the task of super-resolution of turbulent flows in the Rayleigh-Benard convection problem. Across a diverse set of evaluation metrics, we show that MeshfreeFlowNet significantly outperforms existing baselines. Furthermore, we provide a large scale implementation of MeshfreeFlowNet and show that it efficiently scales across large clusters, achieving 96.80% scaling efficiency on up to 128 GPUs and a training time of less than 4 minutes.Comment: Supplementary Video: https://youtu.be/mjqwPch9gDo. Accepted to SC2

The stochastic counterpart of conservation laws with heterogeneous conductivity fields: application to deterministic problems and uncertainty quantification

Conservation laws in the form of elliptic and parabolic partial differential equations (PDEs) are fundamental to the modeling of many problems such as heat transfer and flow in porous media. Many of such PDEs are stochastic due to the presence of uncertainty in the conductivity field. Based on the relation between stochastic diffusion processes and PDEs, Monte Carlo (MC) methods are available to solve these PDEs. These methods are especially relevant for cases where we are interested in the solution in a small subset of the domain. The existing MC methods based on the stochastic formulation require restrictively small time steps for high variance conductivity fields. Moreover, in many applications the conductivity is piecewise constant and the existing methods are not readily applicable in these cases. Here we provide an algorithm to solve one-dimensional elliptic problems that bypasses these two limitations. The methodology is demonstrated using problems governed by deterministic and stochastic PDEs. It is shown that the method provides an efficient alternative to compute the statistical moments of the solution to a stochastic PDE at any point in the domain. A variance reduction scheme is proposed for applying the method for efficient mean calculations

MeshfreeFlowNet: A Physics-Constrained Deep Continuous Space-Time Super-Resolution Framework

We propose MeshfreeFlowNet, a novel deep learning-based super-resolution framework to generate continuous (grid-free) spatio-temporal solutions from the low-resolution inputs. While being computationally efficient, MeshfreeFlowNet accurately recovers the fine-scale quantities of interest. MeshfreeFlowNet allows for: (i) the output to be sampled at all spatio-temporal resolutions, (ii) a set of Partial Differential Equation (PDE) constraints to be imposed, and (iii) training on fixed-size inputs on arbitrarily sized spatio-temporal domains owing to its fully convolutional encoder. We empirically study the performance of MeshfreeFlowNet on the task of super-resolution of turbulent flows in the Rayleigh-Benard convection problem. Across a diverse set of evaluation metrics, we show that MeshfreeFlowNet significantly outperforms existing baselines. Furthermore, we provide a large scale implementation of MeshfreeFlowNet and show that it efficiently scales across large clusters, achieving 96.80% scaling efficiency on up to 128 GPUs and a training time of less than 4 minutes. We provide an open-source implementation of our method that supports arbitrary combinations of PDE constraints